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The origin of a new kinematically identified metal-poor stellar stream, the KFR08 stream, has not been established. We present stellar 
parameters, stellar ages, and detailed elemental abundances for Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Ni, Zn, Sr, Y, Zr, Ba, La, and Eu for 
16 KFR08 stream members based on analysis of high resolution spectra. Based on the abundance ratios of 14 elements, we use the 
chemical tagging method to identify the stars which have the same chemical composition, and thus, might have a common birthplace, 
such as a cluster. Although three stars were tagged with similar elemental abundances ratios, we find that, statistically, it is not certain 
that they originate from a dissolved star cluster. This conclusion is consistent with the large dispersion of [Fe/H] (cr[Fc/H] = 0.29) 
among the 16 stream members. We find that our stars are a enhanced and that the abundance patterns of the stream members are well 
matched to the thick disk. In addition, most of the stream stars have estimated stellar ages larger than 11 Gyr. These results, together 
with the hot kinematics of the stream stars, suggest that the KFR08 stream is originated from the thick disk population which was 
perturbed by a massive merger in the early universe. 
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1. Introduction 


5—I Star streams, which are groups of stars on the same orbit in the 
c/) Galactic potential, have been detected as over-densities in the 
Cd velocity distribution of stars in the solar neighbourhood (e.g., 
““'Arifyanto & Fuchs 2006 Dehnen l998[ Famaey et al!l2005| l. 


Helmi et al.|(|1999]r discovered the signature of a cold stream in 


ical properties of the Hercules stream are coupled to dynamical 
interactions with the Galactic bar. 

Recently, a new candidate stream, called the KFR08 stream, 
on a quite radial orbit was discovered by Klement et al. (2008i 
exploring the RAVE DRl experimental data ( |Steinmetz et al.| 
|2006|l. This stream is present as a broad feature in the range 


the velocity distribution of the halo and interpreted this stream 
^—1 as part of the tidal debris of a disrupted satellite galaxy accreted 
OO by the Milky Way. Moreover, several studies ( |Arifyanto & Fuchs] 
|20061|Helmi et al.|2006||Navarro et al. 12004)1 have concluded that 
OO the Arcturus stream is another such debris stream dating back to 
an accretion event 5 to 8 Gyr ago. The external origin of such 
'sO halo streams is supported by numerical simulations ( [Helmi et al.| 
[2003 1 and observations of ongoing satellites accretion such as 
in that of the Sagittarius dwarf galaxy ( jlbata et al.|1994| ). 

^ However, accreted satellites are not the only source of 
streams. A stellar stream that is homogenous in age and chemi- 
cal composition is associated with a dissolved star cluster. One 
^ striking example is the HR 1614 str eam (jPe Silva et al.|2007bt 


-180 km s ^ < y < -140 km s ' centered at V 


-160 km s 


and would belong to the stellar halo population. The velocities, 
including U and W, here are given relative to the local standard 
of rest (LSR). Because of the high W-velocities of the stars in 
the stream, an origin external to the Milky Way’s disk was sug¬ 
gested b}jKlemenF£ra^^200^. This is supported by [Klementj 
et al.|(20()9|l and Bobylev et al. (2010|l analyzing two indepen¬ 


dent dat sets. [Klement et al.| ( |2009| ) characterized the orbits of 
calibrated stars from SDSS DR7 data ( [Abazajian et al.[[2009| l 
through angular momentum, eccentricity, and orbital polar an¬ 
gle. An overdensity region that corresponds to the location of the 
KFR08 stream within parameter distributions was found. Based 
on a new version of the Hipparcos catalogue ([van Leeuwen 


^ [Eggen 1978[ Feltzing & Holmberg[[2000| l. Stellar streams can 


also originate from dynamical effects within the disks due to res¬ 
onances with the bars or spiral arms ( [Antoja et al.|2009] Dehnen| 
2000 [l. Analysis of high-resolution spectra of nearby F and G 


2007 1 and updated Geneva-Copenhagen Survey (Holmberg et al. 


2007 1 of F and G dwarfs, [Bobylev et al.[ (2010i identified sta 


dwarf stars has revealed that the stars in the Hercules stream 
have a wide range of stellar ages, metallicities and element abun¬ 
dances ([Bensby et al.|2007)l. They concluded that the kinemat- 


* Based on observations made with Nordic Optical Telescope at La 
Palma under programme 44-014 and on data obtained from the ESO 
Science Archive Facility under programme ID 07LB-0529(A), 072.C- 
0488(E), 077.C-0192(A), 082.B-0610(A), 085.C-0062(A). 


tistically significant signals of the main inhomogeneities in the 
velocity distribution using wavelet transform technique. They 
found that 19 stars cluster around (-160, 225) km s * in the V 
and (U^ + 2y^)F2 plane. However, the KFR08 stream could not 
be completely confirmed by re-analyzing the pure RAVE DR2 
dwarf sample ( [Klement et al.[[2011| l. Five stars, at the phase- 
space position of the KFR08 stream, identified from DR2 have a 
large scatter in metallicity. This is inconsistent with the predic¬ 
tion of chemical homogeneity of stream members if they come 
from a single cluster. But it might agree with the idea proposed 
by [Minchev et al.[ ( [2009| l that this high-velocity stream has a dy- 
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namical origin in the thick disc. It is thought to arise due to a 
sudden energy kick imposed by a massive satellite in the past. 
The signature of this perturbation can be identified in the stel¬ 
lar kinematics if it was not wiped out by the radial mixing (see 
|Bland-Hawthorn et al.|2010| l within the thick disc. 

Although brilliant works on dynamical analysis have been 
done, the detailed elemental abundance and age structure of 
the KFR08 stream have never been investigated. Analysis of 
high-resolution spectra of possible stream members could give 
us more clues for exploring the origin of the stream. In this 
work, we present how candidates of the stream were selected 
from databases and how the spectra of them were observed us¬ 
ing different instruments in Section 2. Stellar parameters, ages 
and elemental abundances are determined from the spectra in 
Section 3. In Section 4, we use chemical tagging to identify the 
putative cluster stars from kinematically detected KFR08 stream. 
The possible origins of the stream are discussed in Section 5. Fi¬ 
nally, conclusions are drawn in Section 6. 


2. Observations 

2. 1. The stellar sample 

[Bobylev et al.| ( |20T0| l found nineteen stars that are likely mem¬ 
bers of the KFR08 stream. Of those 17 are high probably mem¬ 
bers. For 16 of these we have analyzed high resolution spectra. 
These stars and their properties are listed in Table We also 
give their Galactic coordinates, distances to the Sun and space 
velocities with respect to the LSR in Table 

2.2. Spectroscopic observations 

FIES observation: Observations were carried out at the Nordic 
Optical Telescope (NOT) for 10 of the candidates using the fibre- 
fed Echelle Spectrograph (FIES) on January 10 - 12 in 2012. A 
solar spectrum was also obtained by observing the sky at day¬ 
time. The wavelength range of the spectra is 370 - 730 nm, with 
a resolving power R ~ 67000 and the average signal-to-noise 
ratio (S/N) is larger than 100 per pixel for most of the spec¬ 
tra. All the spectra were reduced using the EIEStooQpipeline. 
The pipeline includes the following steps to reduce the observed 
frame: subtracting bias and scattered light, dividing by a nor¬ 
malized two-dimensional flat field, extracting individual orders, 
and solving wavelength. Finally, all individual spectral orders 
are merged into a one-dimensional spectrum. 

Archival data: Reduced ID spectra for 6 stars were extracted 
from the ESO archive. We found that 3 stars have been observed 
using the UVES spectrograph ( |Dekker et al.|2000| l on the VET 
8-m telescope in 2013. Two spectra have medium resolution R 
~ 45000, while the spectrum of star HIP 59785 was obtained in 
high resolution mode (R ~ 110000). The S/N values of those 
three spectra are larger than 100 per rebined pixel. The spec¬ 
trum for HIP 117702 was observed in 2006 with EEROS on ESO 
2.2m-telescope ( |Kaufer et al.|1999) l. This spectrum has a large 
wavelength range (350-920 nm), high resolution (R ~ 48000) 
and good S/N (~ 75 per rebined pixel). The exoplanet survey 
carried out with HARPS on nearby stars with high resolution 
spectroscopy (R ~ 120000, [Mayor et al.||2003] l includes spec¬ 
tra for two of the stars identifled as possible stream members in 
[Bobylev et al. ( [2010[ l. Both spectra have a mean S/N larger than 
50 per rebined pixel. 

* http://www.not.iac.es/instruments/fies/fiestool/FIEStool.html 
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Radial velocities (RVs) were measured by cross-correlating 
the solar synthesis spectrum with the observed spectra, making 
use of the packages RVSAO and XCSAO within IRAp[^ Both our 
own spectra and collected spectra from the archive are shifted to 
rest wavelength by cross-correlation with a synthetic solar spec¬ 
trum using the IRAE task DOPCO. 


3. Abundances Analysis 

3.1. Stellar Parameters 


Eor our spectral analysis, we use Spectroscopy Made Easy 
(SME, [Valenti & Fisdier|[2005i [Valenti & Piskunov[[T996l l to 
determine both stellar parameters and elemental abundances by 
comparing synthetic spectra with observations. In SME the 
model atmospheres are interpolated in the precomputed MARCS 
model atmosphere grid ( [Gustafsson et al. 20081, which have 
standard composition. 

The stellar parameters (Teff, togg, [Ee/H]) were determined 
following the same methodology as in Liu et al.| ( 2015[ hereafter 
L15). We refer the reader to LI5 for full details. In brief, initial 
stellar parameters which are used to generate a synthetic spec¬ 
trum were estimated before running SME. As most of the can¬ 
didates have metallicities around -0.7 dex (Table 2 in [Bobylev[ 
[etal.|2010| l, we assume an initial [Ee/H] of -0.7 dex for all stars. 
Eollowing the same methods as in LI5, initial effective tem¬ 
peratures (Tgff) were determined using the colour-metallicity- 
temperature relations of Alonso et al. ( 1996| l, using both B - V 
and b-y colors, where B-V came from the Hipparcos catalogue 
([Perryman et al.|1997[l and b - y values were obtained from the 


Stromgren survey ( 01sen|1983 |1993[[Sdiuster & Nissen|1988] l. 

An initial estimate of surface gravities (log g) were calculated us¬ 
ing Teff, distances (from Hipparcos parallaxes), bolometric cor¬ 
rections, and stellar mass from the Yonsei-Yale isochrones ( [De-[ 
[marque et al.|2004] l. Then, the stellar parameters (shown in Ta- 
ble[ljl were determined by using a purely spectroscopic method 
(called Procedure 1 in L15). The location of stars in Teff vs. 
log g diagram are shown in Pig[T] Considering systematic errors 
and possible sources of uncertainty in the atmospheric model 
and atomic line data, the systematic uncertainties in the stel¬ 
lar parameters were estimated to be 6T^s - 67 + 40 K, dlog 
g = 0.08 + 0.06 dex, and d[Pe/H] = 0.06 + 0.03 dex. 

Since our stellar parameters are derived using Ee lines, which 
were compiled for the spectra of solar type stars, it is possible 
that the results may include systematic uncertainties due to de¬ 
partures from non-local thermodynamical equilibrium (NLTE, 
e.g.|Bergemann et al.|2012|l. Eollowing the same methodology 


as in [Ruchti et al. ( 2013[ ), we derived stellar parameters for two 
stars (HIP 5336 and HIP 58708) using on the fly NLTE correc¬ 
tions for Ee. It was found that the mean differences between our 
results and stellar parameters determined from equivalent width 
method are ATeff = -67 K, Alog g = 0.02 dex and A [Ee/H] = 
-0.02 dex, respectively. These differences are well within our 
estimated uncertainties, and will not effect our chemical tagging 
experiment. 

As the precision of the parallax for all stars is better than 
15%, based on the partly physical method (called Procedure 2 
in L15), we also measured their stellar parameters (Table [T]) by 
fixing logg as estimated from the Hipparcos parallaxes. Com¬ 
parison with the logg estimated from the Ee lines, the differ- 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy (AURA) under cooperative agreement with the National Sci¬ 
ence Fundation. 
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ences between the two surface gravities are less than 0.15 dex 
for most of the stars. An exception is found for two stars, HIP 
59785 and HIP 87101, their loggr are changed by 0.6 dex. It 
was also found that the mean differences in Teu are larger than 
200 K for those two stars. Since HIP 59785 is a very cool star 
with B — V > 1.0, we might overestimate log g from parallax 
caused by the wrong bolometric correction. Because this cool 
star is close to the colour limits for the bolometric correction 
(see Torres|[2010[ ). For HIP 87101, it is apparent that the loca¬ 
tion of this star is far from the isochrones in the Fig. How¬ 
ever, the Teff and log g based on the parallaxes are coupled with 
the isochrones. It should also be notice that HIP 87101 is the 
most metal-poor star in our sample. It might be subject to strong 
NLTE effects (see |Bergemann et al.|2012| ). These imply that we 
might obtain um'eliable stellar parameters from the Procedure 1 
for this star. We noticed that the difference in [Fe/H] is -0.28 
dex for the star HIP 55988. The noised Fe lines might induce 
this huge difference in [Fe/H] because of the lowest S/N (< 30 
per rebined pixel) spectrum among the stars. Except for these 
three stars, the mean difference between two obtained stellar pa¬ 
rameters are ATeff = 59 + 66 K, Alogp = 0.12 + 0.16 dex, and 
A[Fe/H] = 0.03 + 0.04 dex, respectively. Since a little changes 
in [Ee/H] and Teff do not affect our final results from the chem¬ 
ical tagging experiments (see Sect. 4). We can conclud that two 
previous derived stellar parameters are consistent at certain level. 
Eor the remainder of our analysis we adopt the stellar parameters 
derived using a purely spectroscopic method. 



Fig. 1. Teff vs. log g diagram. Filled circles indicate our main sam¬ 
ple from Hipparcos catalogue. Isochrones at different ages (10 and 13 
Gyr) and metallicities (-0.6 and -1.0 dex) according to the Yonsei-Yale 
models ([Demarque et al.|2004[) are shown. 


3.2. Stellar Ages 

Using the stellar parameters derived from the spectra, stellar ages 
were estimated from fits to the Yonsei-Yale isochrones ( |Demar-| 
que et al.|2004| l by maximising the probability distribution func¬ 
tions as described in [Bensby et al.| ( |201 1[ ). The most probable 
age is determined from the peak of the age probability distribu¬ 
tion, Icr lower and upper age limits are obtained from the shape 
of the distribution. Stellar masses were also determined in a sim¬ 


ilar manner. Both the ages and masses are listed in Table It 
is clearly seen that most of the stars, for which ages can be de¬ 
termined, have ages larger than 11 Gyr. The possible system¬ 
atic biases mainly caused by sampling the isochrone data points 
( [Nordstrom et al.|2004| l in our probabilistic age determinations 
was discussed in L15. Comparing with the given typical errors 
in ages, they found that it can be ignored. 


3.3. Elemental Abundances 


Abundances of cr-elements (Mg, Si, Ca, Ti), iron peak elements 
(Cr, Ni, Zn), odd-Z elements (Na, Al), s-process elements (Sr, 
Y, Zr, Ba, La), r-process (Eu), and Scandium were measured by 
fitting the selected absorption lines for each element simultane¬ 
ously. Except for Ba, La, and Eu, we combined the line lists 
from the |Bensby et al. ( |2003| l and the Gaia-ESO line list (Heiter 
et al. in prep) as we did in L15. In this work, the average abun¬ 
dance of Ti and Cr were derived by measuring both neutral and 
singly ionized lines. Since hyperfine splitting (HES) has the ef¬ 
fect of desaturating strong lines (|McWilliam et al.|1995|l, we an 


alyzed the four Ba II lines by adopting the HFS from Me William 
(|1998|l. And for the La and Eu lines, we adopted the data from 
|Lawler et al.| ( |200Ta| l and |Lawler et al.| ( |2001b| l, respectively. In 
addition, isotopes further contribute to the broadening of Ba, La, 
and Eu lines (e.g. [Lawler et al.|2001bHSneden et al.|2002| l. We 
assumed solar system isotopic ratios for these elements. 

During the abundance analysis, we left the corresponding el¬ 
emental abundance (e.g. [Na/H]) free while the stellar parame¬ 
ters were kept fixed. The average ratio ANa (the absolute abun¬ 
dance relative to the total number density of atoms) is the output 
from SME. The solar elemental abundance values taken from 
[Grevesse et al.| ( [2007[ l were used to normalize the SME results, 
in order to obtain the abundance ratio, e.g., [Na/H]. The abun¬ 
dance ratios with respect to Fe (i.e. [X/Fe] in standard notation) 
were also calculated and are shown in Fig. and and in Ta¬ 
ble To normalize the determined abundances of our sample, 
we also derived solar abundances, using the same line list, and 
the stellar parameters derived from our solar spectrum in Sect. 
3.1. 


SME gives a typical error in abundance less than 0.01 dex 
due to continuum placement and line blending. However, the 
uncertainties in our derived abundances are dominated by the un¬ 
certainty in stellar parameters. In Sect. 3.1 the uncertainties in 
the stellar parameters were found to be o-t^^ - 40 K, triogg - 0.06 
dex, and cr[pe/H] = 0.03 dex. The uncertainties in the elemental 
abundances associated with these, for two stars (HIP 18235, and 
HIP 64920), were calculated. The total uncertainty, shown in Ta¬ 
ble]^ was derived by taking the square root of the quadratic sum 
of the different errors, as L15 suggested. The average values of 
the total uncertainties for all elements are between 0.03 and 0.05 
dex. Although more than two lines for the r-and s-process ele¬ 
ments were carefully selected from the literature, there are some 
cases that only a single line was used to derive the elemental 
abundance for a given star. This could introduce more uncer¬ 
tainty to the abundance than what is expected. 


4. Chemical tagging 

The chemical tagging method described by [Mitschang et al. 
( 2013[ l was employed to find out whether the KFR08 stream 
originates from a dissolved stellar cluster. Here we give a brief 
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Table 1. Stellar parameters of KFR08 stream members. 


Names 

HIP 

Tcff 

( K ) 

log 9 

[ Fe / H ] 

^mic 

( kms -*) 

^mac 

( kms -') 

vsini 

( kms -‘) 

T' 

-'eff 

( K ) 

logg' 

[ Fe / H ]' 

v' ■ 

mic 

(km s “') 

^mac 

(km s “’) 

vsini' 
(km s “^) 

5336 

5304 

4.52 

- 0.77 

1.0 

2.0 

0.7 

5354 

4.60 

- 0.76 

0.7 

2.0 

0.6 

15495 

5029 

4.60 

- 0.57 

0.6 

1.0 

0.6 

5119 

4.79 

- 0.55 

0.7 

0.7 

0.8 

18235 

4980 

3.33 

- 0.68 

0.6 

3.2 

1.0 

4976 

3.30 

- 0.71 

1.1 

3.3 

0.7 

19143 

4803 

4.76 

- 0.88 

1.0 

0.9 

0.5 

4939 

5.07 

- 0.84 

0.5 

0.3 

0.4 

54469 

5948 

4.20 

- 0.45 

1.0 

3.0 

2.7 

6093 

4.42 

- 0.42 

1.2 

1.8 

3.5 

55988 

4096 

4.65 

- 0.82 

0.4 

- 0.1 

0.1 

4071 

5.08 

- 0.54 

0.2 

0.0 

0.1 

58357 

4898 

3.20 

- 0.67 

0.6 

3.3 

0.9 

4912 

3.25 

- 0.68 

1.1 

3.4 

0.8 

58708 

4869 

3.41 

- 0.60 

1.0 

3.0 

0.9 

4924 

3.37 

- 0.56 

1.1 

3.1 

0.8 

58843 

5649 

4.37 

- 0.76 

1.0 

2.9 

0.5 

5636 

4.39 

- 0.79 

1.0 

2.9 

0.5 

59785 

4620 

1.98 

- 0.63 

1.5 

4.9 

1.9 

4808 

2.52 

- 0.48 

1.7 

4.7 

1.5 

60747 

5268 

4.54 

- 0.85 

0.5 

1.5 

0.6 

5319 

4.65 

- 0.84 

0.6 

1.6 

0.6 

64920 

4672 

2.47 

- 0.58 

1.2 

0.9 

4.1 

4701 

2.51 

- 0.57 

1.2 

0.7 

4.3 

74033 

5657 

4.03 

- 0.73 

1.1 

3.5 

1.3 

5646 

4.02 

- 0.75 

1.1 

3.3 

1.8 

81170 

5069 

4.23 

- 1.22 

0.8 

0.0 

0.1 

5280 

4.78 

- 1.13 

0.5 

- 0.2 

0.1 

87101 

5485 

3.14 

- 1.68 

1.5 

- 0.2 

0.1 

5875 

4.23 

- 1.43 

1.3 

- 0.5 

0.1 

117720 

5141 

4.58 

- 0.56 

0.6 

0.5 

1.5 

5150 

4.63 

- 0.65 

0.8 

1.0 

0.8 


Notes. The second to seventh columns give the global parameters obtained through purely spectroscopic approach (Procedure 1, LI5), while the 
eighth to thirteenth columns give the six parameters obtained from Procedure 2 in LI5. The unprimed quantities are the preferred values which 
are used to estimate the stellar masses and ages. 


Table 2. The Galactic coordinates, distances, space velocities, radial velocities, masses, absolute magnitudes and ages of KFR08 stream members. 


Names 

HIP 

1 

( deg ) 

b 

( deg ) 

d 

( pc ) 

o-d 

( pc ) 

U " 

( kms -‘) 

V “ 

(km s^L 

W " 
(km s “*) 

Vr 

(km s “*) 

O'Vr 

(km s “^) 

M 

1 Mo 

My 

( mag ) 

O' My 

( mag ) 

Age 

( Gyr ) 

-IfT 

( Gyr ) 

+ 1{T 

( Gyr ) 

5336 

125.3637 

- 07.8683 

7.6 

0.1 

-32 

-153 

-28 

- 96.2 

0.2 

0.70 

5.76 

0.01 

14.7 

9.9 

— 

15495 

155.1775 

- 19.8313 

46.2 

2.8 

58 

-174 

-3 

- 104.5 

0.2 

0.69 

6.34 

0.13 

13.1 

6.3 

13.5 

18235 

217.8451 

^ 8.7113 

64.2 

2.4 

-16 

-161 

-19 

120.6 

0.5 

0.87 

2.65 

0.08 

14.7 

9.3 

— 

19143 

163.4802 

- 14.1808 

36.6 

2.3 

-140 

-143 

-42 

109.9 

0.3 

0.61 

7.15 

0.14 

6.4 

3.2 

12.1 

54469 

284.1956 

14.8416 

107.9 

15.0 

91 

-159 

-64 

156.4 

0.5 

0.94 

4.66 

0.30 

8.9 

8.1 

10.5 

55988 

254.0288 

62.2288 

27.4 

1.3 

50 

-154 

-25 

23.6 

0.4 

0.50 

8.02 

0.10 

14.6 

7.8 

— 

58357 

145.6480 

66.4998 

146.0 

21.1 

-123 

-134 

45 

47.2 

0.1 

0.87 

2.53 

0.31 

14.7 

9.9 

— 

58708 

151.6559 

71.2086 

57.1 

1.7 

-14 

-160 

15 

- 10.6 

0.1 

0.87 

2.84 

0.06 

14.8 

11.5 

— 

58843 

275.4489 

63.6586 

70.2 

6.6 

122 

-138 

-58 

10.3 

0.1 

0.77 

4.99 

0.20 

14.8 

12.7 

— 

59785 

295.8082 

20.4609 

123.2 

6.8 

-117 

-136 

-109 

25.2 

1.0 

1.92 

0.8 

0.12 

0.6 

— 

6.5 

60747 

288.9989 

63.7756 

81.9 

10.7 

no 

-146 

91 

153.4 

0.2 

0.68 

5.91 

0.28 

14.6 

9.7 

— 

64920 

310.1019 

34.1637 

165.6 

15.6 

66 

-159 

43 

140.3 

0.5 

0.88 

0.70 

0.20 

11.0 

5.7 

13.1 

74033 

9.8705 

53.0007 

71.9 

6.3 

-113 

-132 

42 

- 59.4 

0.2 

0.83 

3.98 

0.19 

14.5 

11.9 

— 

81170 

11.6638 

27.7129 

45.1 

2.8 

-77 

-157 

-123 

- 170.9 

0.1 


6.34 

0.13 




87101 

16.9365 

9.4726 

110.0 

15.4 

-76 

-159 

-3 

- 129.5 

0.3 

1.22 

4.48 

0.30 

1.9 

1.8 

7.0 

117720 

315.0151 

- 54.2925 

50.4 

2.7 

12 

-159 

124 

- 25.4 

0.9 

0.70 

5.93 

0.12 

13.1 

6.7 

13.6 


Notes. indicates the Icr lower (or upper) age could not he determined because a very young (or old) star is out side of the isochrones limitations. 
“ means that the space velocities of the stars come from|Bobylev et aL](|2010^. 


Table 3. The elemental abundances of KFR08 stream members. 


HIP 

Na 

Mg 

Al 

Si 

Ca 

Sc 

Ti 

[ X / Fe ] 

Cr 

Ni 

Zn 

Sr 

Y 

Zr 

Ba 

La 

Eu 

5336 

0.05 

0.22 

0.27 

0.16 

0.25 

0.10 

0.31 

0.04 

- 0.02 

0.11 

- 0.03 

- 0.09 

0.22 

- 0.10 

0.01 

0.29 

15495 

0.25 

0.20 

0.42 

0.14 

0.35 

0.09 

0.49 

0.13 

0.00 

0.03 

- 0.02 

- 0.01 

0.30 

- 0.16 

- 0.08 

0.40 

18235 

0.09 

0.36 

0.31 

0.23 

0.28 

0.10 

0.34 

0.04 

- 0.01 

0.12 

0.13 

0.11 

0.37 

0.08 

0.09 

0.37 

19143 

0.27 

0.38 

0.47 

0.06 

0.49 

0.14 

0.70 

0.21 

- 0.01 

0.00 

0.06 

0.01 

0.75 

- 0.02 

0.31 

0.51 

54469 

0.11 

0.29 

0.06 

0.17 

0.16 

0.10 

0.22 

0.02 

0.03 

0.20 

0.36 

0.33 

0.36 

0.14 

0.02 

0.21 

55988 

0.36 

0.30 

0.49 

0.43 

0.25 

0.14 

0.34 

0.03 

0.01 

0.29 

- 0.19 

0.01 

- 0.17 

- 0.23 

0.39 

0.68 

58357 

0.12 

0.33 

0.35 

0.25 

0.26 

0.13 

0.35 

0.07 

0.00 

0.10 

0.29 

0.27 

0.43 

0.45 

0.46 

0.35 

58708 

0.21 

0.33 

0.44 

0.33 

0.28 

0.19 

0.37 

0.10 

0.08 

0.14 

0.08 

0.03 

0.24 

- 0.10 

- 0.01 

0.35 

58843 

0.03 

0.35 

0.20 

0.20 

0.21 

0.06 

0.25 

0.00 

- 0.04 

0.09 

0.06 

0.00 

0.20 

0.00 

0.05 

0.35 

59785 

0.20 

0.43 

0.43 

0.45 

0.21 

0.17 

0.19 

0.02 

0.01 

0.16 

0.08 

- 0.10 

0.04 

- 0.32 

- 0.19 

0.37 

60747 

0.01 

0.31 

0.29 

0.15 

0.24 

0.02 

0.31 

0.09 

- 0.05 

0.05 

0.02 

0.04 

0.27 

0.08 

- 0.07 

0.47 

64920 

0.21 

0.36 

0.38 

0.29 

0.23 

0.16 

0.31 

0.03 

0.05 

0.11 

0.20 

- 0.05 

0.06 

- 0.10 

- 0.15 

0.24 

74033 

0.09 

0.36 

0.20 

0.18 

0.21 

0.12 

0.27 

0.00 

- 0.03 

0.11 

0.12 

0.03 

0.34 

0.10 

0.10 

0.26 

81170 

- 0.16 

0.30 

0.13 

0.13 

0.41 

- 0.18 

0.35 

0.05 

- 0.09 

0.14 


0.12 


- 0.12 


- 0.12 

87101 

- 0.07 

0.53 

- 0.27 

0.30 

0.35 

- 0.18 

0.20 

- 0.22 

0.06 

0.08 


- 0.10 

0.05 

- 0.04 

- 0.20 

0.10 

117720 

- 0.04 

0.45 

0.37 

0.14 

0.44 

0.09 

0.51 

0.14 

- 0.04 

0.07 

0.22 

0.15 

0.51 

- 0.11 


0.39 


Notes. Abundances of 16 elements (Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Ni, Zn, Sr, Y, Zr, Ba, La, and Eu) relative to Fe are listed in columns 2 to 17. 
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Fig. 2. Light elemental abundance ratios [X/Fe] as the function of [Fe/H]. Three tagged cluster stars (see Sect. 4) are indicated in red filled 
diamonds, while other stream stars are marked by black emptied diamonds. The disk stars from |Bensby et al.]^2014[l are marked open gray circles, 
while the cyan triangles and magenta stars indicate the thick disk and halo stars collected from |Nissen & Schuster] ( [2010^ and |Fulbright| ( |2000l l, 
respectively. Dashed lines indicate solar values. 


Table 4. Total random errors in the abundances due to the uncertainties in stellar parameters. 



_ 

_ 

_ 

_ 

_ 

_ 

_ 

o-[X/Fe] 

_ 

_ 

_ 

_ 

_ 

_ 

_ 

_ 

Name 

Na 

Mg 

Al 

Si 

Ca 

Sc 

Ti 

Cr 

Ni 

Zn 

Sr 

Y 

Zr 

Ba 

La 

Eu 

HIP 18235 

±0.03 

±0.05 

±0.04 

±0.03 

±0.05 

±0.03 

±0.05 

±0.04 

±0.04 

±0.03 

±0.03 

±0.05 

±0.05 

±0.03 

±0.04 

±0.04 

HIP 64920 

±0.04 

±0.04 

±0.03 

±0.03 

±0.05 

±0.03 

±0.05 

±0.02 

±0.02 

±0.03 

±0.02 

±0.03 

±0.03 

±0.03 

±0.03 

±0.02 


summary of this method. A metric (5c) was defined as 

Sc - / 
c 


Nr 


( 1 ) 


where Nc is the number of measured abundances, A^ and A^ 
are individual abundance ratios of element C with respect to Fe 
relative to solar for stars i and j, respectively. When the element 
C is Fe, then Ac is the ratio of Fe to H. Here, wc is the weight¬ 
ing factor for an individual species. For our purposes this is set 
to unity as we do not know if any of the elements are more or 
less important. Furthermore, 5c is the mean absolute difference 
between any two stars across all measured elements. To turn this 
5c into a probability that shows how likely it is that the stars 
come from a dissolved cluster, we use an empirecal probability 
function from |Mitschang et al.| ( |2013[ ) that turns the 5c into a 
probability Ps^. 

Given a confidence limit P/™, we first removed all pairs for 
which the probability is less than this threshold. Secondly, all the 
remaining candidates that make up the pairs were re-evaluated 
based on their 5c and Pg^. In this way we have separated the 
potential dissolved clusters stars from the field stars using their 
elemental abundances. As Mitschang et al.|p013| l suggested, the 
cluster detection confidence, Pdus, is finally evaluated using the 
mean of 5c in the tagged group. 

Two confidence limits 85% and 68% were used to verify 
whether the stream comes from a dissolved cluster or not. The 


68% threshold, which is analogous to a Icr detection, is the low¬ 
est meaningful probability and corresponds to a 5c of 0.058 dex 
which is comparable with the systematic uncertainty in the abun¬ 
dances. As [Mitschang et ar] ( |2013| l point out, a linked group of 
stars are less contaminated by other groups if a higher confidence 
level is used (P//,„ = 85%). 

We can not obtain abundances of [Sr/Fe], [Zr/Fe] and [La/Fe] 
for three stars, because the absorption lines are weak. Thus abun¬ 
dances of 14 elements (Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Fe, Ni, 
Zn, Y, Ba, and Eu) are used to identify the stream members. We 
found that three stars, HIP 18235, HIP 58843, and HIP 74033 
could belong to one group using a limit of P/„„ = 68%. The 
mean 5c of the three stars is 0.044 which gives a group detec¬ 
tion confidence of ~ 84%. No star was tagged as a member of 
a group of stars when we tested with a higher threshold (P;,„, = 
85%). 

In Fig. 1^ we show the abundance patterns of the three tagged 
cluster stars in comparison with the other stream stars. The fig¬ 
ure shows that the three stars (shown as red filled diamonds) al¬ 
most have the same elemental abundances within the uncertain¬ 
ties. If all stars originated from a single star cluster, we would 
expect that they exhibit similar abundance ratios in all elements. 
This can be seen for some elements, such as Mg, Cr, and Ni. The 
scatter of abundance ratios among our sample stars for those ele¬ 
ments are between 0.04 and 0.08 dex which are comparable with 
or relatively larger than the measurement uncertainty. However, 
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Fig. 3. Heavy elemental abundance ratios [X/Fe] as the function of [Fe/H]. The symbols and colors have the same meaning as the Fig.|^ Dashed 
lines indicate solar values. 


the stars have very a large scatter (crjx/pe] ^ 0.12 dex) in abun¬ 
dance ratios in other elements, such as Na, Al, Ti, Y, Ba, and 
Eu. The metallicity, [Fe/H], exhibited the largest dispersion with 
o‘[Fe/H] = 0.29. These large dispersions suggest that the stars 
have different birthplaces. As shown in Fig. the star-to-star 
scatter of chemical composition among three tagged stars is sig¬ 
nificantly different from the other stars. This illustrates that our 
chemical tagging method, based on a selection of 14 elements, 
is efficient to isolate the cluster stars from the field stars. 

As mentioned before, the stellar parameters based on the par¬ 
allaxes were computed using the Procedure 2 for our sample. 
Then, we measured their elemental abundances using derived 
stellar parameters in the same way as introduced in Sect. 3. The 
same three stars were identified as cluster members when we 
followed the same chemical tagging experiment as before. 


5. Discussion 


5.1. Is the KFR08 a dissolved cluster? 


Recent studies on observations (see |De Silva et al.|2007a Pan-| 
|cino et ar]|2010| ) and simulations (Feng & Krumholz 2014) l 
demonstrated that chemical compositions within star clusters are 
homogenous. The stream stars should have the same elemen¬ 
tal abundances and stellar age if the KFR08 stream originated 
from a cluster. As mentioned in Sect. 4, the star-to-star scat¬ 
ter in abundances of most elements are quite large compared to 
the measurement uncertainty. The sample stars spread in abun¬ 
dances space rather than cluster together. We also found that 4 
out of 16 stars (see Table are much younger than the other 
members. These suggest that at least a fraction of the stream 
members are not from a dissolved cluster. 

How likely is it that our three stars found as belonging to 
a chemically homogenous cluster in Sect. 4 really come from a 
dissolved cluster and are not a chance grouping? [Mitschang et al.| 
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( |2013| l predicts that half of the members in a chemically tagged 
group in fact are interlopers when a confidence limit of 68% is 
used. Thus, it is possible that the three stars found in Sect. 4 do 
not belong to a chemically homogenous group of stars, i.e., a dis¬ 
solved cluster. This is supported by our higher confidence level 
experiment that no star is linking to a group when the threshold 
is set to 85%. 

In addition, if we assumed that the three stars, HIP 18235, 
HIP 58843, and HIP 74033, are siblings of a dissolved cluster, 
they should have similar RVs. However, we found that the differ¬ 
ence in their RVs is larger than 50 km s ' (see Table]^. We also 
might expect the three stars to have similar kinematics if they 
were born in the cluster. When we looked at their angular mo¬ 
menta (see Fig.|^ and total velocities (see Fig.|^, however, they 
are not more clustered than other stars in the angular momentum 
and velocity spaces. These further weaken their likelihood of 
being a dissolved stellar cluster. 

The probability that the star HIP 74033 belongs to the 
KFR08 stream is 0.65 according to Bobylev et al. ( 2010| l. Since 
the location of this star in velocity space is also close to the Arc- 
turus stream, we are not surprised that it was selected as a poten¬ 


tial member of the Arcturus stream in Arifyanto & Fuchs (2006)1. 


Thus the nature of this star and which stream it might belong to 
is rather unclear. 

Although we found that three stream candidates have simi¬ 
lar elemental abundances, from this discussion we conclude that 
they do not belong to a dissolved cluster. 


5.2. Does the KFR08 originate from an accreted galaxy? 

A possible origin for the KFR08 stream is that it is part of 
the tidal debris from a satellite accreted long ago. Such a mi¬ 
nor merger origin was proposed by Klement et al. ( |2008) l. A 
significant population of stars, which lag behind the LSR by 
~ 100 km s ', with kinematics intermediate between the canon- 
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Fig. 4. Relative elemental abundance ratios with respect to the Sun ([X/Fe]) for all stars as a function of atomic number. When the element is Fe, 
the relative elemental abundance is [Fe/H]. The red dots with red lines represent the abundance patterns of the three stars tagged as cluster stars 
according to our chemical tagging technique (see Sect. 4), while circles with dotted lines represent the abundance patters of the other stars in our 
sample. The dashed line indicates the solar abundance. The gray break in the x-axis is due to the cut of blank space between Y and Ba elements. 


ical thick disk and the stellar halo were interpreted as remnants 
of a merger by Gilmore et al.| ( |2002| l and |Wyse et’aL] ( |2006| l. 
[Gilmore et al.| ( |2002| ) also found more evidence that the stellar 
halo retains kinematic substructure indicative of minor mergers. 
Our sample has V velocity ^ -152 km s ' with respect to the 
LSR. The stream candidates selected by [Klement et al.j ( |2008| 
from the RAVE survey look like halo stars and show quite radial 
orbits (high W velocities). The mean W velocity of our sample 
is W = -7 + 65 km s“'. We calculated the angular momenta of 
our stream members. These are shown in Fig.|^ We found that 
the stars cluster around = 500 kpc km s“'. This suggests that 
they might belong to the halo population. We also find that the 
stars have large scatter in the component. Thus, 

they do not show any difference when we compare them with the 
distribution of normal halo stars in L, and Lj_ space (e.g.,|Kepley| 
jet al.|20071 l. 


Chemically, our stars have [a/Fe] ratios that are enhanced 
for any given a-element. This is distinct from the low ratios 
observed in surviving dwarf satellite galaxies (e.g., [Tolstoy et al.] 
[2009[ also see Sect. 5.3). Thus, if the KFR08 stars were accreted, 
they must have come from a merger with a more massive satellite 
galaxy. In this case, the stars could have smaller W velocities, 
like our observations, since they would be dragged into the plane 
of the Milky Way by dynamical frication (e.g., [Read et al.|2008[ 
[Ruchti et al.|2014) . 


In addition to the dynamical arguments given above, Venn 
[et al.[ ( [2004 [ l demonstrate that the thick disk cannot be com¬ 
prised of the remnants from a low-mass dwarf spheroidal (dSph) 
galaxy, because dSph galaxies are cr-poor when compared to the 
Galactic stars of similar [Fe/H]. As can be seen in Fig. and 
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Fig. 5. Distribution of the stream candidates in angular momentum (L^ 
and Lj^)-space. Note that the distance of the local standard of rest to the 
Galactic center is Rq = 8.5 kpc. The coloured diamonds have the same 
meaning as in Fig. Two dash lines are used to illustrate the angular 
momenta to typical halo stars. 


Fig. the stream candidates are a enhanced and [Fe/H] of the 
stars cover a broad range. 

5.3. Does the KFR08 have a dynamical origin? 

For the s-and r-process elements, such as Ba, La, and Eu, a large 
scatter in abundance ratios is seen in Fig.[^ Since they are pro¬ 
duced in different places and time scales, this implies that our 
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candidates might have different birth places. We collected a large 
sample consisting of thin and thick disk and halo stars from the 
literature ( |Bensby et al.|201^|Fulbright|2000 Nissen & Schus- 


[ter| 2010| l and show these stars together with our data in Fig. 
and As can be seen, the abundance patterns of our samp’ 
stars are well matched to the thick disk. It has been shown that 
the Eu abundance follows the Mg abundance in thin and thick 
disk stars while halo stars show an overabundance of Eu relative 


to Mg ( Mashonkina & Gehren|[200T Mashonkina et al.||2003| l. 
Eor our sample stars, the mean value [Eu/Mg] = -0.02 + 0.20 is 
consistent with the trends found for the thick disk stars. 

Generally, stars with a total velocity Dtot = {U^ + V^ + 
greater than ~ 70 kms^', but less than ~ 200 km s“\ are 
likely to be thick disk stars (e.g., |Nissen|[2004|l. According to 
data presented in the Toomre diagram (Eigure |^, the kinemat¬ 
ics of the KEROS stream is intermediate between the thick disk 
and halo populations. Combining with the chemical signature of 
thick disk for the stream stars, they might have once belonged 
to the thick disk and gained hot kinematics as a result of a satel¬ 
lite merger. This is consistent with the hypothesis suggested by 
Minchev et al. ( 2009|l that the high-velocity KEROS stream has a 


dynamical origin presumably due to a strong perturbation in the 
Galactic disk from a merger. 

As mentioned before, our sample stars have a large disper¬ 
sion in W velocity. The U velocities also have a large scatter 
with (Tu ~ S6 km s"'. We have shown that the KEROS stream 
is older than 11 Gyr. Such old populations with high velocity 
dispersion could be expected if they were perturbed by massive 
mergers in the early universe ( [Minchev et al.|20T4) l. It is possi¬ 
ble that the stars were also subsequently migrated from the inner 
disk to their current position. However, it is unclear how much 
radial migration would increase the velocity dispersion of the 
stars ( [Minchev et al.|20T3]|Vera-Ciro et al.|2014| l. 
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Fig. 6. Toomre diagram. The two circles delineate constant total 
space velocities of n,o, = 70 and 200 km s“*, respectively. The coloured 
diamonds have the same meaning as in Fig.[^ 


6. Conclusions 


We derived the stellar parameters and elemental abundances for 
16 of the KEROS stream members identified by Bobylev et al. 
( [2010| l by comparing synthetic spectra with observed spectra. 
The stellar ages are also obtained by fitting to isochrones. To 
find out whether the stream is a dissolved cluster, a chemical 
tagging method was used to tag the cluster members. We found 
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that three stars have similar abundances and could belong to one 
group using a confidence limit of = 68%. However, further 
tests (see Sect. 5.1) do not support the conclusion that the three 
chemically similar stars come from a dissolved cluster. 

Except for a dissolved cluster origin, the stream could orig¬ 
inate from an accreted satellite galaxy or have a dynamical ori¬ 
gin caused by a massive merger. Klement et al.[ ( |2008] l propose 
that the KEROS stream is part of the tidal debris from a satel¬ 
lite accreted long ago. Although the stream stars cluster around 
angular momentum - 500 kpc km s ', they have a large scat¬ 
ter in the component. This make them similar to the field 
halo stars in L, and space. The mean U and W velocities of 
our sample are very small with respect to the LSR. This might 
be inconsistent with the hypothesis that the KEROS stream was 
accreted from a dwarf satellite galaxy. 

We further found that the stream members are a-enhanced. 
This speaks strongly against the accretion debris origin of stream 
if we expect the KEROS stream progenitor to be similar to cur¬ 
rent day dSph galaxies. It has been noted that a high-mass dwarf 
galaxy has higher rate of star formation than the remaining satel¬ 
lite of the Milky Way and their metallicity can be enriched to 
[Ee/H] = -0.6 ( [Tamura et al.[[20(JT] ). We thus can not rule out 
a satellite progenitor of the stream which has a substantial mass 
similar to the EMC or Sgr dwarf galaxy. We found that the abun¬ 
dance patterns of stream stars are well matched to the thick disk, 
especially for the [Eu/Mg]. In addition to the very old ages (> 11 
Gyr), the stream stars have hotter kinematics than the canonical 
thick disk stars. A more likely scenario is thus that the members 
of the KEROS stream have a dynamical origin due to a strong 
perturbation from a merger event in early universe. 
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